Introduction to the snow vole dataset

For these exercises we will use some published data, freely available from the Dryad repository: Bonnet, T & Postma, E (2018). Data from: Fluctuating selection and its (elusive) evolutionary consequences in a wild rodent population. Dryad Digital Repository. https://datadryad.org//resource/doi:10.5061/dryad.6767m.

For some context, these data are from the long-term monitoring of a population of snow voles (Chionomys nivalis) in the Swiss Alps.

A young snow vole. Picture by T. Bonnet

A young snow vole. Picture by T. Bonnet

We will use only the file “Phenotypic data” or “YearPheno.txt”.

What exactly the variables are is not crucial, but will help make sense of the exercises:

  • ID: focal individual unique identifier
  • animal: duplicate of ID to fit relatedness matrices in animal model
  • Mother: mother of the focal individual
  • Year: year of measurement
  • Sex: Female or male
  • Age: Juvenile or adult
  • BMI: Body mass index. Check the publication for details
  • RJst: Standardized relative Julian day
  • RJ2st: Squared Standardized relative Julian day
  • Phi and PhiZ: survival to the next year (duplicated for historical reasons)
  • RhoZ: reproduction on the next year
  • FitnessZ: RhoZ + 2Phi. Check the publication for details
  • GGImm: proportion of assumed immigrant ancestry used to fit genetic groups. Check the publication for details
  • BMIst: BMI standardized within years
  • BMI1 and BMI2: BMI in years of positive or negative selection.

Understanding the data

Linear models

Generalized linear models

Mixed models

Generalized linear mixed models

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OgogIGh0bWxfZG9jdW1lbnQ6CiAgICBkZl9wcmludDogcGFnZWQKICAgIHRvYzogeWVzCiAgICB0b2NfZmxvYXQ6IHllcwogIGh0bWxfbm90ZWJvb2s6CiAgICB0b2M6IHllcwogICAgdG9jX2Zsb2F0OiB5ZXMKLS0tCgojIyBJbnRyb2R1Y3Rpb24gdG8gdGhlIHNub3cgdm9sZSBkYXRhc2V0CgpGb3IgdGhlc2UgZXhlcmNpc2VzIHdlIHdpbGwgdXNlIHNvbWUgcHVibGlzaGVkIGRhdGEsIGZyZWVseSBhdmFpbGFibGUgZnJvbSB0aGUgRHJ5YWQgcmVwb3NpdG9yeToKQm9ubmV0LCBUICYgUG9zdG1hLCBFICgyMDE4KS4gRGF0YSBmcm9tOiBGbHVjdHVhdGluZyBzZWxlY3Rpb24gYW5kIGl0cyAoZWx1c2l2ZSkgZXZvbHV0aW9uYXJ5IGNvbnNlcXVlbmNlcyBpbiBhIHdpbGQgcm9kZW50IHBvcHVsYXRpb24uIERyeWFkIERpZ2l0YWwgUmVwb3NpdG9yeS4gW2h0dHBzOi8vZGF0YWRyeWFkLm9yZy8vcmVzb3VyY2UvZG9pOjEwLjUwNjEvZHJ5YWQuNjc2N21dKGh0dHBzOi8vZGF0YWRyeWFkLm9yZy8vcmVzb3VyY2UvZG9pOjEwLjUwNjEvZHJ5YWQuNjc2N20pLgoKRm9yIHNvbWUgY29udGV4dCwgdGhlc2UgZGF0YSBhcmUgZnJvbSB0aGUgbG9uZy10ZXJtIG1vbml0b3Jpbmcgb2YgYSBwb3B1bGF0aW9uIG9mIHNub3cgdm9sZXMgKF9DaGlvbm9teXMgbml2YWxpc18pIGluIHRoZSBTd2lzcyBBbHBzLgoKIVtBIHlvdW5nIHNub3cgdm9sZS4gUGljdHVyZSBieSBULiBCb25uZXRdKGJyYWluLmpwZykKCgoKV2Ugd2lsbCB1c2Ugb25seSB0aGUgZmlsZSAiUGhlbm90eXBpYyBkYXRhIiBvciAiWWVhclBoZW5vLnR4dCIuCgpXaGF0IGV4YWN0bHkgdGhlIHZhcmlhYmxlcyBhcmUgaXMgbm90IGNydWNpYWwsIGJ1dCB3aWxsIGhlbHAgbWFrZSBzZW5zZSBvZiB0aGUgZXhlcmNpc2VzOgoKKiBJRDogZm9jYWwgaW5kaXZpZHVhbCB1bmlxdWUgaWRlbnRpZmllcgoqIGFuaW1hbDogZHVwbGljYXRlIG9mIElEIHRvIGZpdCByZWxhdGVkbmVzcyBtYXRyaWNlcyBpbiBhbmltYWwgbW9kZWwKKiBNb3RoZXI6IG1vdGhlciBvZiB0aGUgZm9jYWwgaW5kaXZpZHVhbAoqIFllYXI6IHllYXIgb2YgbWVhc3VyZW1lbnQKKiBTZXg6IEZlbWFsZSBvciBtYWxlCiogQWdlOiBKdXZlbmlsZSBvciBhZHVsdAoqIEJNSTogQm9keSBtYXNzIGluZGV4LiBDaGVjayB0aGUgcHVibGljYXRpb24gZm9yIGRldGFpbHMKKiBSSnN0OiBTdGFuZGFyZGl6ZWQgcmVsYXRpdmUgSnVsaWFuIGRheSAKKiBSSjJzdDogU3F1YXJlZCBTdGFuZGFyZGl6ZWQgcmVsYXRpdmUgSnVsaWFuIGRheQoqIFBoaSBhbmQgUGhpWjogc3Vydml2YWwgdG8gdGhlIG5leHQgeWVhciAoZHVwbGljYXRlZCBmb3IgaGlzdG9yaWNhbCByZWFzb25zKQoqIFJob1o6IHJlcHJvZHVjdGlvbiBvbiB0aGUgbmV4dCB5ZWFyCiogRml0bmVzc1o6IFJob1ogKyAyUGhpLiBDaGVjayB0aGUgcHVibGljYXRpb24gZm9yIGRldGFpbHMKKiBHR0ltbTogcHJvcG9ydGlvbiBvZiBhc3N1bWVkIGltbWlncmFudCBhbmNlc3RyeSB1c2VkIHRvIGZpdCBnZW5ldGljIGdyb3Vwcy4gQ2hlY2sgdGhlIHB1YmxpY2F0aW9uIGZvciBkZXRhaWxzCiogQk1Jc3Q6IEJNSSBzdGFuZGFyZGl6ZWQgd2l0aGluIHllYXJzCiogQk1JMSBhbmQgQk1JMjogQk1JIGluIHllYXJzIG9mIHBvc2l0aXZlIG9yIG5lZ2F0aXZlIHNlbGVjdGlvbi4KCiMjIFVuZGVyc3RhbmRpbmcgdGhlIGRhdGEKCiMjIExpbmVhciBtb2RlbHMKCiMjIEdlbmVyYWxpemVkIGxpbmVhciBtb2RlbHMKCiMjIE1peGVkIG1vZGVscwoKIyMgR2VuZXJhbGl6ZWQgbGluZWFyIG1peGVkIG1vZGVscwo=